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Interactions between particles in quantum many-body systems can lead to collective behavior 
described by hydrodynamics. One such system is the electron-hole plasma in graphene near the 
charge neutrality point which can form a strongly coupled Dirac fluid. This charge neutral plasma 
of quasi-relativistic fermions is expected to exhibit a substantial enhancement of the thermal con¬ 
ductivity, due to decoupling of charge and heat currents within hydrodynamics. Employing high 
sensitivity Johnson noise thermometry, we report the breakdown of the Wiedemann-Franz law in 
graphene, with a thermal conductivity an order of magnitude larger than the value predicted by 
Fermi liquid theory. This result is a signature of the Dirac fluid, and constitutes direct evidence of 
collective motion in a quantum electronic fluid. 


Understanding the dynamics of many interacting parti¬ 
cles is a formidable task in physics, complicated by many 
coupled degrees of freedom. For electronic transport in 
matter, strong interactions can lead to a breakdown of 
the Fermi liquid (FL) paradigm of coherent quasiparti¬ 
cles scattering off of impurities. In such situations, the 
complex microscopic dynamics can be coarse-grained to 
a hydrodynamic description of momentum, energy, and 
charge transport on long length and time scales [1]. Hy¬ 
drodynamics has been successfully applied to a diverse 
array of interacting quantum systems, from high mobility 
electrons in conductors [2], to cold atoms [3] and quark- 
gluon plasmas [4] . As has been argued for strongly inter¬ 
acting massless Dirac fermions in graphene at the charge- 
neutrality point (CNP) [5-8], hydrodynamic effects are 
expected to greatly modify transport coefficients as com¬ 
pared to their FL counterparts. 

Many-body physics in graphene is interesting due to 
electron-hole symmetry and a linear dispersion relation 
at the CNP [9, 10]. In particular, the Fermi surface van¬ 
ishes, leading to ineffective screening [11] and the forma¬ 
tion of a strongly-interacting quasi-relativistic electron- 
hole plasma, known as a Dirac fluid [12]. The Dirac fluid 
shares many features with quantum critical systems [13]: 
most importantly, the electron-electron scattering time is 
fast [14-17], and well suited to a hydrodynamic descrip¬ 
tion. A number of exotic properties have been predicted 
including nearly perfect (inviscid) flow [18] and a diverg¬ 
ing thermal conductivity resulting in the breakdown of 
the Wiedemann-Franz law [5, 6]. 

Away from the CNP, graphene has a sharp Fermi sur¬ 
face and the standard Fermi liquid (FL) phenomenology 
holds. By tuning the chemical potential, we may mea¬ 
sure thermal and electrical conductivity in both the Dirac 


fluid (DF) and the FL in the same sample. In a FL, 
the relaxation of heat and charge currents is closely re¬ 
lated as they are carried by the same quasiparticles. The 
Wiedemann-Franz (WF) law [19] states that the elec¬ 
tronic contribution to a metal’s thermal conductivity Kg 
is proportional to its electrical conductivity a and tem¬ 
perature T, such that the Lorenz ratio C satisfies 


C = 


Ke 

ctT 
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( 1 ) 


where e is the electron charge, fee is the Boltzmann con¬ 
stant, and Cq is the Sommerfeld value derived from FL 
theory. Cq depends only on fundamental constants, and 
not on specific details of the system such as carrier den¬ 
sity or effective mass. As a robust prediction of FL the¬ 
ory, the WF law has been verified in numerous metals 

[19] . However, in recent years, an increasing number of 
non-trivial violations of the WF law have been reported 
in strongly interacting systems such as Luttinger liquids 

[20] , metallic ferromagnets [21], heavy fermion metals 
[22], and underdoped cuprates [23], all related to the 
emergence of non-Fermi liquid behavior. 

The WF law is expected to be violated at the CNP 
in a DF due to the strong Coulomb interactions between 
thermally excited charge carriers. An electric field drives 
electrons and holes in opposite directions; collisions be¬ 
tween them introduce a frictional dissipation, resulting 
in a finite conductivity even in the absence of disorder 
[24]. In contrast, a temperature gradient causes electrons 
and holes to move in the same direction inducing an en¬ 
ergy current, which grows unimpeded by inter-particle 
collisions (Fig. 3C inset). The thermal conductivity is 
therefore limited only by the rate at which momentum is 
relaxed due to residual impurities. 
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FIG. 1. Temperature and density dependent electrical and thermal conductivity. (A) Resistance versus gate voltage 
at various temperatures. (B) Electrical conductivity (blue) as a function of the charge density set by the back gate for different 
bath temperatures. The residual carrier density at the neutrality point (green) is estimated by the intersection of the minimum 
conductivity with a linear fit to log((T) away from neutrality (dashed grey lines). Curves have been offset vertically such that 
the minimum density (green) aligns with the temperature axis to the right. Solid black lines correspond to 4e^/h. At low 
temperature, the minimum density is limited by disorder (charge puddles). However, above Tdis ~ 40 K, a crossover marked 
in the half-tone background, thermal excitations begin to dominate and the sample enters the non-degenerate regime near 
the neutrality point. (C-D) Thermal conductivity (red points) as a function of (C) gate voltage and (D) bath temperature 
compared to the Wiedemann-Franz law, cjTCo (blue lines). At low temperature and/or high doping (|/r| 3> feeT), we find the 
WF law to hold. This is a non-trivial check on the quality of our measurement. In the non-degenerate regime (|p| < IcbT) 
the thermal conductivity is euhanced and the WF law is violated. Above Tei-ph ~ 80 K, electron-phonon coupling becomes 
appreciable and begins to dominate thermal transport at all measured gate voltages. All data from this figure is taken from 
sample S2 (inset lE). 


Realization of the Dirac fluid in graphene requires that 
the thermal energy be larger than the local chemical po¬ 
tential /r(r), defined at position r: /cbT > lAi(r)|. Impu¬ 
rities cause spatial variations in the local chemical po¬ 
tential, and even when the sample is globally neutral, it 
is locally doped to form electron-hole puddles with finite 
/i(r) [25-28]. Formation of the DF is further complicated 
by phonon scattering at high temperature which can re¬ 
lax momentum by creating additional inelastic scattering 
channels. This high temperature limit occurs when the 
electron-phonon scattering rate becomes comparable to 
the electron-electron scattering rate. These two temper¬ 
atures set the experimental window in which the DF and 
the breakdown of the WF law can be observed. 

To minimize disorder, the monolayer graphene samples 
used in this report are encapsulated in hexagonal boron 
nitride (hBN) [29]. All devices used in this study are 
two-terminal to keep a well-defined temperature profile 


[30] with contacts fabricated using the one-dimensional 
edge technique [31] in order to minimize contact resis¬ 
tance. We employ a back gate voltage Vg applied to 
the silicon substrate to tune the charge carrier density 
n = rie — u-h, where rie and nu are the electron and hole 
density, respectively (see supplementary materials (SM)). 
All measurements are performed in a cryostat controlling 
the temperature Tbath- Fig. lA shows the resistance R 
versus Vg measured at various fixed temperatures for a 
representative device (see SM for all samples). From this, 
we estimate the electrical conductivity a (Fig. IB) using 
the known sample dimensions. At the CNF, the residual 
charge carrier density rimin can be estimated by extrap¬ 
olating a linear fit of log((T) as a function of log(n) out 
to the minimum conductivity [32]. At the lowest tem¬ 
peratures we find rimin saturates to ^8x10® cm“^. We 
note that the extraction of rimm by this method overesti¬ 
mates the charge puddle energy, consistent with previous 
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reports [29]. Above the disorder energy scale T^is "^40 K, 
’T-min increases as Tbath is raised, suggesting thermal ex¬ 
citations begin to dominate and the sample enters the 
non-degenerate regime near the CNP. 

The electronic thermal conductivity is measured us¬ 
ing high sensitivity Johnson noise thermometry (JNT) 
[30, 33]. We apply a small bias current through the sam¬ 
ple that injects a joule heating power P directly into the 
electronic system, inducing a small temperature differ¬ 
ence AT = Tg —Tbath between the graphene electrons and 
the bath. The electron temperature Tg is monitored inde¬ 
pendent of the lattice temperature through the Johnson 
noise power emitted at 100 MHz with a 20 MHz band¬ 
width defined by an LC matching network. We designed 
our JNT to be operated over a wide temperature range 
3-300 K [33]. With a precision of ^ 10 mK, we mea¬ 
sure small deviations of Tg from Tbath, i-S. AT Tbath- 
In this limit, the temperature of the graphene lattice is 
well thermalized to the bath [30] and our JNT setup al¬ 
lows us to sensitively measure the electronic cooling path¬ 
ways in graphene. At low enough temperatures, elec¬ 
tron and lattice interactions are weak [33, 34], and most 
of the Joule heat generated in graphene escapes via di¬ 
rect diffusion to the contacts (SM). As temperature in¬ 
creases, electron-phonon scattering becomes appreciable 
and thermal transport becomes limited by the electron- 
phonon coupling strength [34-36]. The onset temper¬ 
ature of appreciable electron-phonon scattering, Tgi_ph, 
depends on the sample disorder and device geometry: 
Tgi_ph ^80 K [33, 34, 37, 38] for our samples. Below this 
temperature, the electronic contribution of the thermal 
conductivity can be obtained from P and AT using the 
device dimensions (SM). 

Fig. 1C plots Ke(bg) alongside the simultaneously mea¬ 
sured cr(lg) at various fixed bath temperatures. Here, for 
a direct quantitative comparison based on the WF law, 
we plot the scaled electrical conductivity as aTCo in the 
same units as Kg. At low temperatures, T < This 40 K, 
where the puddle induced density fluctuations dominates, 
we find Kg « crTTo, monotonically increasing as a func¬ 
tion of carrier density with a minimum at the neutrality 
point, confirming the WF law in the disordered regime. 
As T increases (T > T^is), however, the measured Kg be¬ 
gins to deviate from the FL theory. We note that this 
violation of the WF law only appears close to the CNP, 
with the measured thermal conductivity maximized at 
n = 0 (Fig ID). The deviation is the largest at 75 K, 
where Kg is over an order of magnitude larger than the 
value expected for a FL. This non-FL behavior quickly 
disappears as jnj increases; Kg returns to the FL value 
and restores the WF law. In fact, away from the CNP, 
the WF law holds for a wide temperature range, consis¬ 
tent with previous reports [33, 34, 37] (Fig. IE). For this 
FL regime, we verify the WF law up to Tei_ph ~ 80 K. 
Finally, in the high temperature regime T > Tgi_ph, the 
additional electron-phonon cooling pathway causes the 
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FIG. 2. Breakdown of the Wiedemann-Pranz law in 
the Dirac fluid regime. The Lorenz ratio is shown as 
a function of the charge carrier density and bath tempera¬ 
ture. Near the CNP and for temperatures above the disor¬ 
der (charge puddle) regime but below the onset of electron- 
phonon coupling, the Lorenz ratio is measured to be an order 
of magnitude greater than the Fermi liquid value of 1 (blue). 
The WF law is observed to hold outside of the Dirac fluid 
regime. All data from this figure is taken from sample SI. 

measured thermal conductivity to be larger than Kg. We 
find that near the CNP Kg tends to decrease just before 
Tei_ph, restricting the maximal observable violation of 
the WF law. 

Our observation of the breakdown of the WF law in 
graphene is consistent with the emergence of the DF. 
Fig 2 shows the full density and temperature depen¬ 
dence of the experimentally measured Lorenz ratio in 
order to highlight the presence of the DF. The blue col¬ 
ored region denotes £ ~ To, suggesting the carriers in 
graphene exhibit FL behavior. The WF law is violated 
in the DF (yellow-red) with a peak Lorenz ratio 22 times 
larger than £o- The green dotted line shows the corre¬ 
sponding rimin (T) for this sample; the DF is found within 
this regime, indicating the coexistence of thermally pop¬ 
ulated electrons and holes. We find that disorder and 
phonon scattering bound the temperature range of the 
Dirac fluid, Tdis < T < Tgi_ph. 

We investigate the effect of impurities on hydrody¬ 
namic transport by comparing the results obtained from 
samples with varying disorder. Fig. 3A shows rimin 
as a function of temperature for three samples used in 
this study. UyainiT — 0) is estimated as 5, 8, and 
10x10® cm“® in samples SI, S2, and S3, respectively. 
All devices show qualitatively similar Dirac fluid behav¬ 
ior; the largest value of £/£o measured in the Dirac fluid 
regime is 22, 12 and 3 in samples SI, S2, and S3, respec¬ 
tively (Fig 3B). For a direct comparison, we show £(n) 
for all three samples at the same temperature (60 K) in 
Fig 3C. We find that cleaner samples not only have a 
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FIG. 3. Disorder in the Dirac fluid. (A) Minimum car¬ 
rier density as a function of temperature for all three sam¬ 
ples. At low temperature each sample is limited by disorder. 
At high temperature all samples become limited by thermal 
excitations. Dashed lines are a guide to the eye. (B) The 
Lorentz ratio of all three samples as a function of bath tem¬ 
perature. The largest WF violation is seen in the cleanest 
sample. (C) The gate dependence of the Lorentz ratio is well 
fit to hydrodynamic theory of Ref. [5, 6]. Fits of all three 
samples are shown at 60 K. All samples return to the Fermi 
liquid value (black dashed line) at high density. Inset shows 
the fitted enthalpy density as a function of temperature and 
the theoretical value in clean graphene (black dashed line). 
Schematic inset illustrates the difference between heat and 
charge current in the neutral Dirac plasma. 


more pronounced peak but also a narrower density de¬ 
pendence, as predicted [5, 6 ]. 

More quantitative analysis of C{n) in our experiment 
can be done by employing a quasi-relativistic hydrody¬ 
namic theory of the DF incorporating the effects of weak 
impurity scattering [5, 6 , 39]. 


C = 


J^DF 

(1 + (n/norf 


( 2 ) 


where 


CuF — 




and 


2 _ min 


( 3 ) 


Here up is the Fermi velocity in graphene, a^nin is the elec¬ 
trical conductivity at the CNF, H is the fluid enthalpy 
density, and Im is the momentum relaxation length from 


impurities. Two parameters in Eqn. (2) are undeter¬ 
mined for any given sample: Im and H. For simplic¬ 
ity, we assume we are well within the DF limit where 
Im and H are approximately independent of n. We fit 
Eqn. (2) to the experimentally measured £{n) for all 
temperatures and densities in the Dirac fluid regime to 
obtain Im and H for each sample. Fig 3C shows three 
representative fits to Eqn. (2) taken at 60 K. Im is esti¬ 
mated to be 1.5, 0.6, and 0.034 ^m for samples SI, S2, 
and S3, respectively. For the system to be well described 
by hydrodynamics, Im should be long compared to the 
electron-electron scattering length of ^ 0.1 fim expected 
for the Dirac fluid at 60 K [18]. This is consistent with 
the pronounced signatures of hydrodynamics in SI and 
S2, but not in S3, where only a glimpse of the DF appears 
in this more disordered sample. Our analysis also allows 
us to estimate the thermodynamic quantity 'H{T) for the 
DF. The Fig. 3C inset shows the fitted enthalpy density 
as a function of temperature compared to that expected 
in clean graphene (dashed line) [18], excluding renormal¬ 
ization of the Fermi velocity. In the cleanest sample H 
varies from 1.1-2.3 eV/^m^ for Tdis < T < Tei_ph. This 
enthalpy density corresponds to ^ 20 meV or ~ 4 A:bT 
per charge carrier — about a factor of 2 larger than the 
model calculation without disorder [18]. 

To fully incorporate the effects of disorder, a hydrody¬ 
namic theory treating inhomogeneity non-perturbatively 
is necessary [40, 41]. The enthalpy densities reported 
here are larger than the theoretical estimation obtained 
for disorder free graphene, consistent with the picture 
that chemical potential fluctuations prevent the sample 
from reaching the Dirac point. While we find thermal 
conductivity well described by Ref. [5, 6 ], electrical con¬ 
ductivity increases slower than expected away from the 
CNF, a result consistent with hydrodynamic transport in 
a viscous fluid with charge puddles [41]. 

In a hydrodynamic system, the ratio of shear viscos¬ 
ity 77 to entropy density s is an indicator of the strength 
of the interactions between constituent particles. It is 
suggested that the DF can behave as a nearly perfect 
fluid [18]: 77 /s approaches a conjecture by Kovtun-Son- 
Starinets: {r]/s)/{h/kB) ^ l/ 47 r for a strongly inter¬ 
acting system [42]. A non-perturbative hydrodynamic 
framework can be employed to estimate 77 , as we discuss 
elsewhere [41]. A direct measurement of 77 is of great 
interest. 

We have experimentally discovered the breakdown of 
the WF law and provided evidence for the hydrodynamic 
behavior of the Dirac fermions in graphene. This pro¬ 
vides an experimentally realizable Dirac fluid and opens 
the way for future studies of strongly interacting rela¬ 
tivistic many-body systems. Beyond a diverging thermal 
conductivity and an ultra-low viscosity, other peculiar 
phenomena are expected to arise in this plasma. The 
massless nature of the Dirac fermions is expected to re¬ 
sult in a large kinematic viscosity, despite a small shear 
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viscosity 77 . Observable hydrodynamic effects have also 
been predicted to extend into the FL regime [43]. The 
study of magnetotransport in the DF will lead to further 
tests of hydrodynamics [5, 39]. 
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Supplementary 

Materials 

SAMPLE FABRICATION 

Single layer graphene is encapsulated in hexagonal 
boron nitride on an n-doped silicon wafer with 285 nm 
Si 02 [31] and is subsequently annealed in vacuum for 15 
minutes at 350 °G. It is then etched using reactive-ion¬ 
etching (RIE) to define the width of the device. A second 
etch mask is then lithographically defined to overlap with 
the sample edge, leaving the rest of the sample rectangu¬ 
lar shaped with the desired aspect ratio. After the RIE 
is performed, the same etch mask is used as the metal 
deposition mask, upon which Gr/Pd/Au (1.5 nm / 5 nm 
/ 200 nm) is deposited. The resulting Ohmic contacts 
show low contact resistances and small PN junction ef¬ 



SI 

S2 

S3 

length (gm) 

3 

3 

4 

width (pm) 

9 

9 

10.5 

mobility (10® cm^ • • s“^) 

3 

2.5 

0.8 

Wmin (10® cm“®) 

5 

8 

10 


TABLE S.I. Basic properties of our three samples. 


fects due to their minimum overlap with device edge. 


OPTIMIZING SAMPLES FOR HIGH 
FREQUENCY THERMAL CONDUCTIVITY 
MEASUREMENTS 

To measure the electronic thermal conductivity Ke of 
graphene using high frequency Johnson noise the sample 
design should be made with three additional considera¬ 
tions: stray chip capacitance, resistance of the lead wires, 
and sample dimensions that enhance electron diffusion 
cooling over phonon coupling. 

Johnson noise thermometry (JNT) relies on measur¬ 
ing the total noise power emitted in a specified frequency 
band and relating that to the electronic temperature on 
the device; to maximize the sensitivity, high frequency 
and wide bandwidth measurements should be made [33] . 
In the temperature range discussed here, the upper fre¬ 
quency limit for JNT is typically set by the amount of 
stray capacitance from the graphene, lead wires, and con¬ 
tact pads to the Si back gate. This is minimized by us¬ 
ing short, narrow lead wires and small (50 pm x 50 pm) 
bonding pads resulting in an estimated 4 pF stray capac¬ 
itance. 

The amount of Johnson noise emitted between any two 
terminals is proportional to the mean electronic tempera¬ 
ture between them where each point in space is weighted 
by its local resistance. Therefore, to maximize the signal 
coming from the graphene, contact resistance should be 
kept at a minimum. To compensate for the narrow lead 
wires, we deposit a thicker layer (200 nm) of gold result¬ 
ing in an estimated total contact resistance of < 80 Q. 

Lastly, to effectively extract Ke from the total elec¬ 
tronic thermal conductance Gth we want to enhance the 
electron diffusion cooling pathway with respect to the 
electron-phonon cooling pathway (see below). This can 
be accomplished by keeping the length of the sample 
short as the total power coupled into the lattice scales 
as the area of the device while diffusion cooling scales as 
1/i?. In addition, the device should be made wide to min¬ 
imize the effects of disordered edges. We find these high 
aspect ratio samples (~ 3:1) are ideal for our measure¬ 
ments and serve the additional purpose of lowering the 
total sample resistance allowing us to impedance match 
over a wider bandwidth. 
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FIG. SI. 2-terminal resistance R vs. back gate voltage for 
the 3 samples used in this report. 


DEVICE CHARACTERIZATION 



FIG. S2. Reflectance TZ = |5iip for a graphene device 
impedance matched to 50 near 125 MHz. 


we can see two critical components of JNT: impedance 
matching over a wide bandwidth and low noise amplifi¬ 
cation. 


In this study we measure three graphene devices en¬ 
capsulated in hexagonal boron nitride (hBN), whose ba¬ 
sic properties are detailed in Table S.I. All devices are 
two-terminal with mobility estimated as 


L 

neRW" 


(S.I) 


where L and W are the sample length and width re¬ 
spectively, e is the electron charge, and n is the charge 
carrier density. The gate capacitance per unit area 
Cg « 0.11 fF/pm^ is estimated considering the 285 nm 
Si02 and ~ 20 nm hBN dielectrics. From this we esti¬ 
mate the charge density 


- Vd) 
e 


(S.2) 


where Vd is the gate voltage corresponding to the charge 
neutrality point (CNP) estimated by the location of the 
maximum of the curve R{Vg). Fig. SI shows the resis¬ 
tance of all samples as a function of gate voltage. 


JOHNSON NOISE THERMOMETRY 


The full Johnson noise thermometry (JNT) setup used 
in this study is outlined in detail in [33]. Here, we only 
give a brief synopsis for completeness. 

The electronic transport within a dissipative device can 
be determined by the high frequency noise power col¬ 
lected by a low noise amplifier as 


{V^) = fesTe X Re(Z)A/ 


1 - 


( Z-Zq 
\Z + Zo 


2 


(S.3) 


where Z is the complex impedance of the device under 
test, Zq is the impedance of the measure circuit (typi¬ 
cally 50 O) and A/ is the bandwidth. From this formula. 


Graphene devices have a typical channel resistance on 
the order of /i/4e^ ^ 6 kO near the CNP. To compensate 
for this, we use an inductor-capacitor (LC) tank circuit 
mounted directly on the sample package to impedance 
match the graphene to the 50 H measurement network. 
Fig. S2 shows the reflectance coefficient TZ = IS'np for a 
typical graphene device after impedance matching. The 
bandwidth and measurement frequency of our JNT is set 
by the Q-factor and LC time of this matchingnetwork. 

At 10 K, the power emitted by a resistor in a 1 Hz 
bandwidth is ~ 10“^^ W or —190 dBm. To amplify 
this signal we use a SiCe low noise amplifier (Caltech 
CITLF3) with a room temperature noise figure of about 
0.64 dB in our measurement bandwidth, corresponding 
to a noise temperature of about 46 K. We operate the am¬ 
plifier at room temperature, outside of the cryostat, to 
ensure it is unaffected by the 3-300 K temperature ramp 
used for thermal conduction measurements. After am¬ 
plification, a homodyne mixer and low pass filter define 
the measurement bandwidth and the power is found by 
an analog RF multiplier operating up to 2 GHz (Analog 
Devices ADL5931). The result is a voltage proportional 
to the Johnson noise power which - after calibration - 
measures the electron temperature in the graphene de¬ 
vice. 

Calibration of our JNT device must be done on every 
sample as each device has a unique R(T, Vg) and there¬ 
fore couples differently to the amplifier. The graphene 
device being measured is placed on a cold finger in a 
cryostat with varying temperature Tbath- With no exci¬ 
tation current in the graphene, we collect the JNT signal 
Vs{T, Vg) for all temperatures and gate voltage needed in 
the study, as shown in Fig. S3. The linear temperature 
slope at each point gives a gain factor g{T, V^) = dVs/dT. 
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FIG. S3. Output voltage K from the JNT measuring a 
graphene device with no excitation current. This is used to 
calibrate the JNT to a given sample. 


MEASURING ELECTRONIC THERMAL 
CONDUCTANCE 



Heating Current (|aA) 

FIG. S4. The electronic temperature in an encapsulated 
graphene device as a function of heating current for three 
different bath temperatures. Te = Tbath + I^R/Gth- The to¬ 
tal thermal conductance between the electronic system and 
the bath is found through Fourier’s Law (solid lines). 


The procedure used to measure electronic thermal con¬ 
ductance is outlined in [30, 33, 44]. A small sinusoidal 
current /(w) is run through the graphene sample, caus¬ 
ing a Joule heating power P(2w) to be injected directly 
into the electronic system. This causes a small tempera¬ 
ture difference AT{2uj) between the electronic tempera¬ 
ture and the bath, described by Fourier’s law: 

P = GthAT. (S.4) 

Here Gth is the total thermal conductance between the 
electronic system and the bath. The component of John¬ 
son noise at frequency 2uj is measured by a lock-in ampli¬ 
fier and then converted to a temperature difference AT 
using the gain g(T, Vg) described in the previous section. 
Fig. S4 shows Te as a function of heating current I for a 
graphene device at three different bath temperatures: 3, 
30 and 300 K. 

THERMAL MODEL OF GRAPHENE 
ELECTRONS 


( P = |2R ) 

4 - 4 



FIG. S5. Simplified thermal diagram of the electronic cooling 
pathways in graphene relevant for our experimental condi¬ 
tions. A current induces a heating power into the electronic 
system which conducts to the bath via two parallel pathways; 
diffusion and coupling to phonons. 


In the regime presented here, Gth is dominated by two 
electronic cooling pathways. Hot electrons can diffuse 
directly out to the contacts (Gdis), or they can couple to 
phonons (Gei-ph): 

G « Gdiff + Gel-ph- (S-5) 

In a typical metal, electron diffusion is described by the 
WF law which is linear in Te. The electron-phonon cool¬ 
ing pathway has two components: first, the electrons 
must transfer heat to the lattice via electron-phonon cou¬ 
pling, and then the lattice must conduct the heat to the 


bath. Consistent with previous experimental studies on 
graphene [33, 34, 37, 38], we find Gei-ph is bottlenecked 
by the weak electron-phonon coupling and hence the lat¬ 
tice is well thermalized to the bath. Fig. S5 shows the 
simplified thermal diagram of the electronic cooling path¬ 
ways in graphene, relevant to our experiment. 

At low temperature, Gth is dominated by Gdig, while 
at high temperature it is dominated by Gei-ph- Fig. S6 
shows an illustration of this effect in one of our devices. 
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FIG. S6. The total thermal conductance of a graphene device 
in the high density regime, illustrating how Gth is dominated 
by either Gdift ~ T, as predicted by the WF law, or Gei-ph ~ 


MEASURING THERMAL CONDUCTIVITY 


In the diffusion-limited regime, we extract the elec¬ 
tronic thermal conductivity Ke as follows. This is detailed 
in [30] and we review the basic calculation for clarity. The 
total power dissipated, P, is given by 
J2 

P=—LW = PLW (S.6) 

<7 

where L is the length of the sample in the direction of cur¬ 
rent flow, W is the width in the perpendicular direction 
and P is the local power dissipated. Because this calcu¬ 
lation is done in linear response, and the external heat 
baths on either side of the sample are at the same temper¬ 
ature, the contributions to power dissipated ^ k(AT)^ 
do not enter so long as AT ~ is small. This is an 
appropriate assumption in the regime of linear response, 
where J is treated as a perturbatively small parameter. 
Fig. S4 shows our experiment is in this regime. 

Let us now determine the change in the temperature 
profile. For simplicity we assume that the graphene sam¬ 
ple is homogeneous, that the approximately uniform elec¬ 
trical current is given by 


, dU dT 

J = -a- - a — , 

ax ax 

and that the heat current is given by 

Q = -aT- -Ke3—, 

ax ax 


(S.7) 


(S.8) 


where 


Ta^ 

Kg = Ke-\ -= Ke(l + ZT). (S.9) 



FIG. S7. Gartoon illustrating the non-uniform temperature 
profile within the graphene-hBN stack during Joule heating 
in the diffusion-limited regime. 


dT /dx is the temperature gradient in the sample, and 
—dV/dx is the electric field in the sample, a/cr is the 
Seebeck coefficient. We assume that the response of 
graphene is dominated only by the changes in voltage 
V and temperature T to a uniform current density J, 
which is applied externally. We also assume that devia¬ 
tions from constant V and T are small, so that the linear 
response theory is valid. Joule heating leads to the fol¬ 
lowing equations: 


-S' 

(S.lOa) 

cr da; ’ 

(S.lOb) 

which can be combined to obtain 


-r. d^T 

(S.ll) 

assuming that Kg is approximately 

homogeneous 


throughout the sample. 

The contacts in our experiment are held at the same 
temperature T. Thus, writing 

T(x) =Te +AT(x), (S.12) 

we find that 

AT{x) = ^—x{L — x). (S.13) 

2/^0 

The average temperature change in the sample, which is 
directly measured through JNT, is 

(AT) = /!FArw = g. (S.14) 

0 

This non-uniform temperature profile is illustrated in 
Fig. S7. Combining Eqs. (S.4), (S.6) and (S.14) we 
obtain 

127 

(S.15) 


In the latter equation, ZT is the thermoelectric coeffi¬ 
cient of merit. As a » 0 at the CNP, we expect ZT « 0, 
and that Rg ~ Kg. 


As we have pointed out in the main text, our samples 
are not perfectly homogeneous, but have local fluctua¬ 
tions in the charge density. Nevertheless, we do recover 











9 


the WF law in the FL regime, suggesting that our mea¬ 
surement of Gth “ and thus Ke - using JNT, along with 
the above formalism, is valid. 


* pkim@physics.harvard.edu 
^ kc.fong@bbn.com 
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